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Abstract. We present here a calculation of the configuration averaged optical 
conductivity of random binary alloys CuAu and AgAu. Our formulation is based on 
the augmented space formalism proposed by Mookerjee [J. Phys. C : Solid State Phys. 
6 1340 (1973)] and the optical conductivity is obtained directly through a recursive 
procedure suggested by Viswanath and Muller ["The user friendly recursion method", 
Troisieme Cycle de la Physique, en Suisse Romande (1993)]. 

PACS numbers: 71.23.-k 
1. Introduction 

In an earlier work pQ we had developed a methodology for the calculation of configuration 
averaged optical conductivity of a disordered alloy based on the augmented space 
formalsim developed by one of us j2j- The method was based on the Kubo-Greenwood 
formula and we had shown that disorder scattering renormalizes both the Green 
functions and the current terms. 

In linear response theory, at zero temperature, the generalized susceptibility of a 
disordered alloy is given by the Kubo formula : 



J — oo 

where, A u (t) is the vector potential, and 

X ^(t-t') = m) &( t -t') ([mnm 

j M is the current operator and is the heaviside step function. If the underlying lattice 
has cubic symmetry, = X The fluctuation dissipation theorem then relates the 
imaginary part of the Laplace transform of the generalized susceptibility to the Laplace 
transform of a correlation function : 

X» = (l/2h) (1 - exp{-/3hu}) S(u) 
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where, 

S(u) = | o °° dt exp{z( W + i5)t} Tr (j"(t) j"(0) ) (1) 

Since the response function is independent of the direction label \i for cubic 
symmetry, in the following we shall drop this symbol. In case of other symmetries 
we have to generalize our results for different directions. Our goal will be, given a 
quantum "Hamiltonian" H to obtain the correlation function, 

s(t) = mmm) 

We shall now determine the correlation directly via the recursion method as 
described by Viswnath and MiillerjS]. In order to simplify the expressions for the 
dynamical quantities as produced by the Hamiltonian, we consider henceforth the 
modified Hamiltonian H = H — E I, whose band energy is shifted to zero. Let 

w)\ = (<t>\m 

The time evolution of this bra is goverened by the Schrodinger equation 

= W*)|H (2) 

We shall now generate an orthogonal basis {{fk\} for representation of equation (J2J. 
We shall do this in the following way : 

(i) We begin with initial conditions : 

(/-il = ; (/ o | = (0ij(O) 

(ii) We now generate the new basis members by a three term recurrance relationship : 

(/ fc+1 | = (/ fc |H - (f k \a k - (f k -i\(3 2 k k = 0,1,2... 



where, 



(f k \H\fk) p 2 (fk\fk) 

ot k = ,, i , . Pk - 



(fk\fk) " (fk-l\fk-l) 

We now expand the bra {i/j(t)\ in this orthogonal basis : 

oo 

<^(*)l = E (M Dkit) (3) 

k=0 

Substituting equation (JHJ) into equation (J2J) and using orthogonality of the basis, we get : 

- iD k {t) = D k _ x {t) + a k D k (t) + D k+l (t) (4) 

with D-i(t) = and D k (0) = 5kfi. We shall now show that the pair of sequences 
generated by us, namely, {ak} and {f3 k } are enough for us to generate the correlation 
function. We note first that : 

D (t) = (m\fo) = S(t) (5) 
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Let us define the Laplace transform : 

dk(z) = dt exp(-izt) DJt) 
Jo 

Putting this back in equation(jIJ) we get : 

(z-a k )d k (z) - iS kfi = d k -i{z) + d k+x {z) k=0,l,2... (6) 
This set of equations can be solved for do( J clS db continued fraction representation : 

M z ) = j2 ( 7 ) 

z-a l - — — 2 

Z — OL\ 



z — a 2 — . .. 

The structure function, which is the Laplace transform of the correlation function 
can then be obtained from the above : 

S(w) = lim 2 do(u + iS) (8) 
2. Augmented space formulation 

Augmented space method was first proposed by one of us [2] as a feasible technique for 
carrying out configuration averaging in disordered systems. In particular, it has been 
applied extensively for electronic structure of disordered alloys. The readers are referred 
to the article [I] for details. Here we shall discuss the salient features of this method for 
application to configuration averages of correlation functions. 

The problem addressed involves the configuration averaging of a function of many 
independent random variables, e.g. <C /({%}) where R labels a set of lattice points. 
The first step is to associate with each random variable ur an operator M fl such that 
the spectral density of this operator is the probability density of the random variable : 

p(^) = -(l/7r)(T|(^I-M^)- 1 |T) 

For example, if ur takes the values and 1 with probabilities x and y = 1 — x then 
the operator M R is of rank 2 and acts on the configuration space of %, <p R , spanned by 
the eigenstates |0) and |1) of M R . The state | f) = y/x\0) + ^/y\l) and the representation 
of M R in the basis {| T>, | I) = y/y\0) - Vx\l)} is 

x 

fxy 

The augmented space theorem then states that : 

«/(W})»=({0}|f({M^})|{0}) (9) 

where f is an operator of rank 2 N obtained by replacing each hr by M^. This operator 
acts on the configuration space of the set of random variables, $ = fl A state 
in this space is given by a pattern of configurations at each site, i.e. sequences like 
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{TT1 • • •} = {C} which can be uniquely represented by the sequence of sites where we 
have (say) a j configuration. This sequence is called a cardinality sequence. The empty 
cardinality sequence, {0}, is the pattern of |-s everywhere. 

If we go back to equations (2)-(5), we notice that for a disordered binary alloy, 
S(t) = 5[H({nij})]. The augmented space theorem then states that : 

« S (t) »=« (0|j(t)j(O)|0> »= (0® {0}|j(t)j(O)|0® {0}) = S[H({M*})] 

where the augmented space Hamiltonian and the current operators are constructed by 

replacing every random variable tir by the corresponding operator M R . 

The recursion may now be modified step by step in the full augmented space : 

m)\ = (0®{0}|J(*) 

The time evolution of this bra is goverened by the Schrodinger equation 

-ij t {m)\) = mm m 

As before, we shall generate the orthogonal basis {{f k \} for representation of 
equation (fTUj) : 

(i) We begin with initial conditions : 

(/-i| = ; (/o| = (0®{0}|j(O) 

(ii) The new basis members are generated by a three term recurrance relationship : 

(/ fc+1 | = (/ fc |H - (f k \a k - (h-tlPt k = 0,l,2... 

where, 

= (/fc|H|/ fc ) ~ 2 = (fk\fk) 
\fk\fk) " (fk-l\fk-l) 

We now expand the bra {ip{t)\ in this orthogonal basis : 

oo 

M*)®{0}| = E (fk\ D k (t) 

k=0 

Continuing exactly as in the last section we get, 

D (t) = (m ® {0}|/ o ) =« S(t) » (12) 
Taking Laplace transforms and using the three term recurrance, 

do(z) = % —^ 2 (13) 

Pi 

z - a " M 

Z — CKi 

z- a 2 - ... 

The configuration averaged structure function, which is the Laplace transform of 
the averaged correlation function can then be obtained from the above : 

< S(uj) > = lim 2 d (co + iS) (14) 

The equations (11)-(14) will form the basis of our calculation of the configuration 
averaged correlation function. 
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3. The Hamiltonian and current operators in augmented space 

Let us begin with a Hamiltonian for a random binary alloy represented in the basis of a 
tight-binding linear muffin-tin orbitals method (TB-LMTO) F° r the sake of 

notational clarity we have suppressed the angular momentum index L. In what follows, 
we may consider R to be a composite index RL. 

H = E Cr Pr + E E ^R 2 Srr' Ar? T rr/ 

R R R' 

Pr and Tr R > are projection and transfer operators on the space spanned by the TB- 
LMTO basis. Cr and A 1 ^ 2 are the TB-LMTO potential parameters and Srr/ the 
structure matrix. For a binary alloy the former are random and may be given by 



C R = C A n R + C B (1 - n R ) =C B + (C A - C B ) n R 

A 1/2 A 1/2 A 1/2 n \ A 1/2 . / A 1/2 A l/2\ 

a r = a a n R + A B ' (l-n R ) = A B / + (A A ' - A B ' ) n R 

Here the random variables hr takes values and 1 with probabilities x and y, as before. 
The augmented space theorem then builds up the augmented Hamiltonian by replacing 
the binary random variables {ur} by the corresponding operators {M R } described in 
the earlier section. 



C R =« C> I + (y - x){C A - C B ) ? r + y^(C A - C B ) T R 

Af =« A 1/2 » I + (y - x)(Af - Af) P l R + V^i^T ~ ^T) ^ (15) 
where the configuration space operators are : 



P R \{C}) = 5(Re{C}) \{C}) 
T$\{C}) = \{C} ± R) 

The augmented space Hamiltonian is then given by : 

H = E C R ® P* + EE ^R 2 SrR'^r? ® Tun, (16) 

R R R' 

Next we look at the expression for the current operator in the TB-LMTO basis : 
\Xr) = l<M + E h RRj\M 

R' 

The dot refers to derivative with respect to energy. In this basis, the matrix elements 
of the current operator can be written as 



j rr> 



(XR'\F\Xr) 



(!)./« 



RR' 



'RR' 



+ E^ 



(2),M 
RR" 



R" 



lR"R> 



RR" Vrhr, 



(3),M 



R" 



EE^ 

R" R" 



RRin Vrihrii It-R"R 



(4),M 



h. 



(17) 
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where 



The technique for calculating these matrix elements has been described in detailed 
by Hobbs et al [6 . We have also used this technique in our earlier paper Tj and we 
shall use it here as well. The readers are referred to these two papers for details. 

Ideally the next step would be to calculate J AA , J^b-- Jab^ Jba as the current terms 
between two sites when they are occupied by atom pairs AA, BB, AB and BA embedded 
in the disordered medium. A simpler first step would be to obtain these current terms 
from the pure A and B and from the ordered AB alloy. In general, the current operator 
can be written as : 

f = E^)P« + EE nx)T RR , (is) 

R R R> 

where x — R ~ R' • 

In a disordered alloy, the current representations are random : 



•HO) = ^a(O) n R + J£ B (0) (1 - n R ) and 
J"{x) = J Ux) n R n R > + JUx) (1 - n R ) (1 - n R >) 

+ J Ux) n R (1 - n R >) + JfiBA(x) (1 - n R ) n R > 

Using the augmented space theorem, we replace the random variable n R by an 
operator M R in the expression of J M (0), J fl (x) we get the current operator in the 
augmented space : 



Y = ( JM( )) I <g, P R + Jf (0) ® P R + j£(0) T# 
+ (J"( x )) I ® T R>R , + (y-x) J[ l ( X ) K ® Tfl, w - 



where 



■ ) 

R! 



R 



■ R,R' 



+ 



%y { Ji(x) 



5Pb 



->T4 



1 T BiB / 
(19) 



< J"(0)> = sJ^O) + yJg B (0) Jf (0) = (y - x)[J^(0) - j£ fl (0)] 
J 2 At (0) = v ^[J^(0)-J^(0)] 

( J(x)> = * 2 J^(x) + ( + JUx)) + y 2 JUx) 

Jf (x) = a: (JJUx) - J5b(x)) + V ( J^( 

J 2 M (x) 
J£(x) 



x (JUx) - JUx)) + y (JUx) - JUx)) 



jH T n jfj, T n 

J AA ~~ J AB ~~ J ™ A ' 



'BA 



BB 
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This augmented current operator is used to construct the starting state of the 
recursion as described in the last section. 

4. Results and Discussion 

For application of our methodology we have chosen two alloy systems disordered : CuAu 
(50-50) and disordered AgAu (50-50). Our choice was governed by two considerations 
: first, we had studied these alloys in an earlier work with a different approach pQ and 
it would be interesting to compare the two approaches ; and secondly, the fact that 
we have experimental data available for both these alloys and it would be interesting 
to compare our theoretical estimates with experiment. We have begun our study with 
the self-consistent TB-LMTO-ASR calculation on AgAu and CuAu 50-50 alloys. We 
have minimized the energy with respect to the variation in the average lattice constant 
for both the alloys. We have calculated the current terms from the basis functions 
and potential parameters obtained from a TB-LMTO calculation for the ordered 50- 
50 alloys. Ideally we should have embedded AA, BB and AB pairs in the disordered 
medium and obtained the current expressions from such a calculation. However, in our 
earlier work we have shown that effect of randomness in the current terms is small. 
This gives us confidence to proceed with these estimates. 

To check the density of states for the disordered alloys we have used our TB-LMTO- 
ASR programme developed by us [Zj . The density of states has been shown in our earlier 
work 1J. The main features of interest is that the filled ci-states lie about 2 eV below 
the Fermi level for both these alloys. 

In figure ^ we show the continued fraction coefficients for a few different angular 
momentum components. The coefficients do show a tendency of showing oscillatory 
convergence, but within the maximum ten steps that we were able to go, it was not 
possible to accurately locate the converged values. We therefore chose not to use the 
square- root terminator, but opted for the Beer-Pettifor terminator [Hj. 

Experimental data on disordered CuAu (50-50) are available jH] • The authors have 
reported high resolution optical data for both ordered and disordered CuAu for photon 
energies between 1.2 eV and 6 eV. In the figure 121 we also show the experimental data 
for e 2 (uj) in dotted lines and compare it with the theoretical data in figure |2fbold lines). 
Drude like behaviour is clearly seen till photon energies of ~ 2.1 eV in theory and 
around 2 eV in the experiment. Subsequently e 2 (uj) rises again because of transitions 
from the ci-states. Both the theory and experiment show a decrease around 6 eV. The 
experimental curve shows a sharper dip and the subsequent peak is at a lower energy. 

We have experimental data on AgAu (50-50) 10 , whose density of states closely 
resembles CuAu. The inter-band contribution to the imaginary part of the dielectric 
function e 2 {oj) may be obtained from the optical conductivity data, by subtracting 
away the Drude contribution and dividing the result by ui : e' 2 (u)) = (cr(^) — a D {uj)) / uj . 
Below the onset of the inter-band transitions, this quantity vanishes. It shows a shoulder 
around 2.5 eV and then reaches a maximum at around 3.5 eV before decreasing and 




showing a feature around 5.5 eV. 

As in the coherent potential approximation (CPA) estimates of Banhart ^0] there is 
a ~ 0.5eV shift outwards in the theoretical estimate. In the above figure we have shifted 
our theoretical data by 0.5eV downwards for comparison. Banhart attributes this shift 
to the foundations of the electronic theory used : local density approximation which 
cannot reproduce excited states properly. This had been discussed also by Fehrenbach 
and Bross 12J. To this we shall add that it is doubtful whether the energy linearized 
TB-LMTO can accurately reproduce the large energy window of our calculation. It may 
be useful in future to base our calculations on the newly developed NMTO formalism 
[T3] . In addition, the termination procedure of the recursion requires further refinement, 
as the convergence of the continued fraction coefficients is not as rapid as in density 
of states calculations. Banhart notes that the CPA neglects all short-ranged order. 
However, since our calculations in the augmented space recursion does include correlated 
scattering from different sites, we cannot attribute the discrepancy to the CPA. 

Apart from the energy shift, the experimental data are in reasonable agreement 
with our theory. The general shape with a maximum is clearly reproduced. The feature 
shown in the experimental data around 5.5 eV is also seen in our theoretical estimate, 
slightly shifted in energy. The shoulder in the experimental data just above 2 eV is not 
properly reproduced in the theoretical curve. The CPA estimate of Banhart (jlD]) also 
does not seem to show such a prominent shoulder. Is it possible that the experimental 
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alloy had short ranged order ? It is known that short-ranged order can have influence 
on optical properties ([14). The augmented state formalism can handle short-ranged 
order ([El) and we shall study this effect in a future communication. 
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Figure 3. Intraband contribution to e^w) for AgAu alloy theoretical (full lines) and 
experimental (Nielsson ^I]) (points) 
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